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1 Introduction 

The Camassa-Holm equation was introduced in [4] in the form 

u t + jXu x — Ujact + 3uu x = 2u x u xx + uu xxx , jU 6 R, (1) 
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and was intensively explored afterwards (see, for example, Refs. 15 6 7, 17 1). Its superiza- 
tions were also constructed, see IT1 I19I . Since 0} is not an evolution equation, its integrabil- 
ity properties (existence and even definition of Hamiltonian structures, conservation laws, 
etc.) are not standard to establish. 

One of the ways widely used to overcome this difficulty is to introduce a new un- 
known m = u — Uxx and transform Eq. {TJ to the system 

m, = -um x - (2m + n)u x , 
m = u — Uxx, 

which has almost evolutionary form. We stress this "almost", because the second equation 
in (O (that can be considered as a constrain to the first one) disrupts the picture and, at best, 
necessitates to invert the operator 1 — D\. At worst, dealing with Eq. ([2]l as with an evolution 
equation may lead to fallacious results. 

In our approach based on the geometrical framework exposed in Ref. |3], we treat the 
equation at hand as a submanifold in the manifold of infinite jets and consider two natural 
extensions of this equation, cf. with Ref 1 16|. The first one is called the ^-covering and serves 
the role of the tangent bundle. The second extension, ^-covering, is the counterpart to the 
cotangent bundle. The key property of these extensions is that the spaces of their nonlocal (in 
the sense of [ 15 1) symmetries and cosymmetries contain all essential integrability invariants 
of the initial equation. The efficiency of the method was tested for a number of problems 
(see Refs. [12 13, 14 1) and we apply it to the Camassa-Holm equation here. 

In Section [2] we briefly expose the necessary definition and facts. Section [3] contains 
computations for the Camassa-Holm equation in its matrix version (computations and re- 
sults are more compact in this representation), while in Section [4] we reformulate them for 
the original form (TJi and compare later the results obtained for the two alternative presenta- 
tions. Finally, Section|5]contains discussion of the results obtained. Throughout our exposi- 
tion we use a very stimulating conceptual parallel between categories of smooth manifolds 
and differential equations proposed initially by A.M. Vinogradov and in its modern form 
presented in Table[T] 



Table 1 Conceptual parallel between two categories 



Manifolds 



Equations 



Smooth manifold 
Point 

Smooth function 
Vector field 
Differential 1-form 
de Rham complex 
Tangent bundle 
Cotangent bundle 



Infinitely prolonged equation 
Formal solution 
Conservation law 
Higher symmetry 
Cosymmetry 

(n — l)st line of Vinogradov's ' 

^-covering 

i* -covering 



-spectral sequence 



This table is not just a toy dictionary but a quite helpful tool to formulate important 
definitions and results. For example, a bivector on a smooth manifold M may be understood 
as a derivation of the ring C°°(M) with values in C°°(T*M). Translating this statement to the 
language of differential equations we come to the definition of variational bivectors and their 
description as shadows of symmetries in the £*-covering (see Theorem [2] below). Another 
example: any vector field (differential 1-form) on M may be treated as a function on T*M 
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(on TM). Hence, to any symmetry (cosymmetry) there corresponds a conservation law on 
the space of the f -covering (£ -covering). This leads to the notions of nonlocal vectors and 
forms that, in turn, provide a basis to construct weakly nonlocal structures (see Subsec- 
tions l3.4l and l3.5l l. Of course, these parallels are not completely straightforward (in technical 
aspects, especially), but extremely enlightening and fruitful. 

The idea of this paper arose in the discussions one of the authors had with Volodya 
Roubtsov in 2007. We agreed to write two parallel texts on integrability of the Camassa- 
Holm equation that reflect our viewpoints. The reader can now compare our results with the 
ones presented in 1 20 1. 

2 Underlying theory 

We present here a concise exposition of the theoretical background used in the subsequent 
sections, see Refs. |3, 13 .131. 



2.1 Equations, symmetries, etc. 

Let n : E — > M be a fiber bundle and K„ : J°°{k) — > M be the bundle of its infinite jets. To 
simplify our exposition we shall assume that n is a vector bundle. In all applications be- 
low % is the trivial bundle R" ! xl"^l". We consider infinite prolongations of differential 
equations as submanifolds S C J°°{%) and retain the notation K„ for the restriction nj\g. 
Any such a manifold is endowed with the Cartan distribution which spans at every point 
tangent spaces to the graphs of jets. A symmetry of S is a vector field that preserves this 
distribution. The set of symmetries is a Lie algebra over R denoted by &ymS. 

For any equation £ its linearization operator lg\ x — > P is defined, where x is the 
moduli of sections of the pullback tz„(tz) and P is the module of sections of some vector 
bundle over S . Then sym<? can be identified with solutions of the equation 

£g((p) = 0, <pex. (3) 

For two symmetries pi, q>2 £ symd 5 their commutator is denoted by {<p\, 92}- 

Denote by A' h the module of horizontal i'-forms on S and introduce the notation 

£ = Homjr(<2,A£), n = dimM, 
for any module Q. The adjoint to lg operator 

ty.p^k (4) 

arises and solutions of the equation 

tg(w) = 0, yeP, (5) 

are called cosymmetries of S\ the space of cosymmetries is denoted by cosymfo 5 . 

Let df, ; A| — > A^ +1 be the horizontal de Rham differential. A conservation law of the 
equation S is a closed form (0 6 A ; " _1 . To any conservation law there corresponds its gener- 
ating function 8(0 £ cosymtf , where 8 : 1 — > E{ n 1 is the differential in the E\ term of 

1 All the modules below are modules over the ring & of smooth functions on S. 
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Vinogradov's "^-spectral sequence, see 1211 . In the evolutionary case 8 coincides with the 
Euler-Lagrange operator. A conservation law is trivial if its generating function vanishes. 
In particular, di, -exact conservation laws are trivial. 

A vector field on £ is called a % '-field if it lies in the Cartan distribution. A differential 
operator A : P — > Q, P and Q being ^-modules, is called a ^-differential operator if it is 
locally expressed in terms of "rf -fields. For example, lg is a ^-differential operator. 

A "if -differential operator H : P — > x is said to be a variational bivector on £ if 

l g oW = W*ot s . (6) 

This condition means that H takes cosymmetries of £ to symmetries. If £ is an evolu- 
tion Eq. l[6} implies also that H* = — H. A bivector H is a Hamiltonian structur^ on £ 
if [[H, H]] s = 0, where [[•, -]] s is the variational Schouten bracket (see also 1111 ). Two Hamil- 
tonian structures are compatible if [[Hi, H2]] s = 0. 

A ^-differential operator S : x — > P is called a variational 2-form on £ if 

t s aS = S*olg. (7) 

Such operators take symmetries to cosymmetries and in evolutionary case are skew-adjoint. 
They are elements of the term E\' n 1 of Vinogradov's "rf -spectral sequence. A variational 
form is a symplectic structure on the equation £ if it is variationally closed, i.e., SS = 0, 
where 8 : E^' n 1 — > E^' n 1 is the corresponding differential. 

We shall also consider recursion ^-differential operators R : x — > x and R : P — > P 
satisfying the conditions 

tgoR = tfolg, t s oR = $!otg (8) 

for some "(f -differential operators R' : P — > P and R' : k — > k. An operator R satisfies the 
Nijenhuis condition if [[R, R]] n = 0, where [[•,-]]" is the variational Nijenhuis bracket. A re- 
cursion operator R is compatible with a Hamiltonian structure H is the composition R o H is 
a Hamiltonian structure as well. 



2.2 Nonlocal theory 

Let £ and £ be equations and f : £ — > £ be a fiber bundle. Denote by *«? and *«? the Cartan 
distributions on £ and £ , resp. We say that | is a covering if for any d € £ the differen- 
tial dg<§ isomorphically maps onto ^tgy A particular case of coverings (the so-called 
Abelian coverings) is naturally associated with closed horizontal 1 -formal 

By definition, any ^-field X on £ can be uniquely lifted to a "if -field X on £ such 
that d^(X) = X. Consequently, any ^-differential operator A: P — > Q is extended to a c €- 
differential operator 

A : # ®& P -> # ®^ g, 
^ being the algebra of smooth functions on <f . 

It is more appropriate to call these objects Poisson structures, but we follow the tradition accepted in the 
theory of integrable systems. 

3 When dimA/ = 2, Abelian coverings are associated with conservation laws of the equation 6°. 
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A nonlocal <§-(co)symmetry of S is a (co)symmetry of the covering equation & . They 
are solutions of the equations £g(p = and t^\jf = 0, resp. Along with these two equations 
one can consider the equations 

(1) l,<p = 0, (2) t g¥ = 0. (9) 

Their solutions are called E, -shadows of symmetries and cosymmetries, resp. A shadow of 
symmetry is a derivation & — > & that preserves the Cartan distributions. For any two shad- 
ows of symmetries (pi and 92 their commutator {91, 92} can be defined. This commutator 
is a shadow in a new covering that is canonically determined by 91 and 92. 

2.3 The I- and £*-coverings 

Let S C J°°{n) be an equation. Its l-covering x: J£(£) — > S is obtained by adding to <f 
the equation £g{q) = 0, where q is a new variable. Dually, the £* -covering x* : J£* {$) — > if 
is constructed by adding the equation = with a new variable p. They are the exact 

counterparts of the tangent and cotangent bundles in the category of differential equations. 
By the reasons that will become clear later, we regard both q and p as odd variables. The 
main point of our method is the fundamental relation between integrability invariants of S 
and shadows in x and T*. To formulate this relation, let us give an auxiliary definition: for 
an arbitrary operator equation Ao A = VoBwe say that A is a trivial solutions if A is of the 
form A = A' oB. Classes of solutions modulo trivial ones will be called nontrivial solutions. 
Then the following results hold. 

Theorem 1 (shadows in the ^-covering) There is a one-to-one correspondence between 
nontrivial solutions of the equation 

igoR = R'oi s , R.x^x, 

and x -shadows of symmetries linear w.r.t. the variables q. In a similar way, there is a one- 
to-one correspondence between nontrivial solutions of the equation 

t g oS = S'o£g, S.x^P, 

and x -shadows of cosymmetries linear w.r.t. the variables q. 

Theorem 2 (shadows in the £* -covering) There is a one-to-one correspondence between 
nontrivial solutions of the equation 

£goH = H'o£*g, H:P->x, 

and X* -shadows of symmetries linear w.r.t. the variables p. In a similar way, there is a 
one-to-one correspondence between nontrivial solutions of the equation 

£* s oR = R! o£* g , R:P^p, 
and T* -shadows of cosymmetries linear w.r.t. the variables p. 

Theorem 3 Let R\ and R2 be recursion operators for symmetries on S and <J"r denote 
the x-shadow corresponding to R. Then \R\ , R2]] 11 = iff {<Pr 1 , ^r, } = 0. Similarly, ifH\ 
and H2 are bivectors then [[Hi,H2]] s =0 iff {^"^ , ^ } =0, where <P^ denotes the X*- 
shadow corresponding to H. 
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In both cases the curly brackets denote the super bracket of shadows that arises due to 
oddness of the variables q and p. Additional discussion of Theorem[3]the reader will find in 
Remark[2] 

Theorem 4 To any cosymmetry of £ there canonically corresponds a conservation law 
of J£{$), Dually, to any symmetry of £ there canonically corresponds a conservation law 

of •/ ■ .: .. 



2.4 Computational scheme 

Let locally the equation £ be given by the system 

(10) 

where j = l,...,m and | a \ <k. 

Step 1 consists of writing out defining equations for symmetries and cosymmetries of £ . 
Let {u J a yj^ s be internal coordinates on £, S and J being some sets of (multi)indices and u{, 
corresponding to d^u-i /dx a . Then any "if -field on £ is a linear combination of the total 
derivatives 

ax aeS.jeJ ou J a 
The linearization of £ is the matrix operator with the entries 

dF l 

(V);= E TT ^ /e'> Z=l,...,r. (12) 

A symmetry <p = (<p , (p m ) enjoys the equation 

dF l 

£ D a (<pi)=0, 1=1,. ..,r, (13) 

aeSJeJ " u a 

and the corresponding field is the evolutionary vector field 

9 9 = £ D a ((pJ)S-, (14) 
aeSJeJ "Kn 



while the bracket of symmetries 91 , 92 is given by 

{<Pb<jP2}' / = 9<s>Mi)-9<p2{^{), j=l,...,m. (15) 
The operator adjoint to ( 112b is 

and a cosymmetry t|/ = (\jf , i/f r ) satisfies the equation 

£(-l) |ff| i3a(f T V , )=0, ye/. (17) 

(7,/ 0«(T 



(^)/=E(-l) |CT| Ocr°— , ye/, / = !,.. .,r, (16) 
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Step 2. Here we look for closed 1 -forms and construct Abelian coverings associated to 
them. A horizontal form CO = Y.iMdx' is closed if 

D x a {A p )=D x p{A a ), 1 < a < jS < n. (18) 

Such a form gives rise to a nonlocal variable w = w m that satisfies the equations 

dw 

= A,-, i=l,...,n. (19) 

ox' 

These equations are compatible on d 1 Ob due to d 1 3 b - Recall that for n = 2 closed 1-forms 
coincide with conservation laws. The total derivatives lifted to the covering equation $ are 

D xi = D x i +Ai^-, i=\,...,n. (20) 
ow 

Step 3. At this step we compute a number of particular symmetries and cosymmetries (using 
equations d 1 3 b and \Y1\ , resp.). They are used to construct canonical nonlocal variables on 
the £*-covering {nonlocal vectors) and on the ^-covering {nonlocal forms), resp., at Step 4. 
We also use them as seed elements in series generated by recursion operators. 



Step 4 consists of construction of the I- and £* -coverings and introduction of canonical 
nonlocal variables over them (see Step 3). The ^-covering is obtained by adding to $9$ the 
system of equations 

„ dF' dqj 
aeSJeJ ou J a dx 
cf. with Eq. d 1 2b , while the £*-covering is given by 

■ dF' I 

dx a 



K- 1 )" 7 ' =°. (22) 

a, I 



that comes from dl6l i. 

If (p is a symmetry of S then one can introduce a covering over ££* (£) described by the 
system 

|f =L A 'oA?)p'o, i=l,...,n, (23) 

a.l 

where A' a ; are ^-differential operators (see Theorem [4}. In a similar way, to any cosymme- 
try y/ there corresponds a covering 

|l=E V iiM«£ i = 1 > ( 24 ) 

V J a j being ^-differential operators as well. We omit here a general description of these 
operators and refer the reader to the particular case of our interest exposed in Sections 14.41 
and 1431 
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Step 5. We now use Theorems [TJand[2]to construct recursion operators and Hamiltonian and 

symplectic structures. Let , . . . , \j/ s cosymmetries of S . Let us consider the covering _£? (&) 
over jSf (<a) with the nonlocal variables q\,...,q s defined by d24b and lift the operators lg 
and tg to this covering. Then the following result specifies Theorem[TJ 

Theorem 5 Let <P = (<P l , <P m ) be a solution of the equation tg(<P) = 0on linear 
w.r.t. q" and qp : 

a.a p 

Then the operator 

R = E a ° JD ° +Jl b p D ^°LKiWDa, 
a p a 

takes shadow of symmetries to shadows of symmetries. In a similar way, to any solution f = 
(f\. ..,W), 

a,c p 

there corresponds the operator 

s = ^c a a jD a +E^' °ZKi(w) D ° 

a p a 

that takes shadows of symmetries to shadows of cosymmetries. 

In a dual way, consider symmetries <pi,...,<p s af the equation S and the covering «Sf * (£*) 
over J£*(£ ) with the nonlocal variables pi,...,p s defined by d23t . Then, lifting lg and t*g, 
we obtain a similar specification of Theorem [2) 

Theorem 6 Let <t> = (& 1 , . . . ,<P m ) be a solution of the equation £g(<P) = on Jzf*(<f) 
linear w.r.t. p„ and pp : 

^=E«" J p"+E^- 

a,a p 

Then the operator 

h = e 4 J Ax + E b i D ~, 1 E 4 ( ) D ° - 

cr y3 cj 

fafes shadow of cosymmetries to shadows of symmetries. In a similar way, to any solu- 
tion W = {W l ,W), 

V J = L c " J P a *+L dj pPP 

a,a p 

there corresponds the operator 

r = E c " + E 4 D * 1 ° £ 4 "•<■ («v» ) D - 

CT J3 CT 

f/iaf ?afej shadows of cosymmetries to shadows of cosymmetries. 

After finding the operators R, S, H and R we check conditions $6^ and ^ and compute 
necessary Schouten and Nijenhuis brackets. 
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Step 6 The last step consists of establishing algebraic relations between the invariants con- 
structed above. 



3 The matrix version 

We consider Eq. {T} in the form 

«r - utxx + 3mm.,' = 2u x u xx + uu xxx , (25) 

i.e., set jii = 0, and, similar to ©, introduce a new variable w = ecu — u xx , where a is a new 
real constant. Consequently, the initial equation transforms to the system 

w t = —2iirW — uw x , 

(26) 

w = ecu — u xx . 

We choose the following variables for internal local coordinates on the infinite prolongation 
of Eq. lESfr : 

d k w d k u d k+l u 

x ,t, Wl = -^, uo* = w ,u 1 t = -^, k = 0, 1, . . . 

Then the total derivatives in these coordinates will be of the form 

D, =3-+ X>*+13— + L M i.*3 + L Z) '(° !M - M; 3 ' 

A =37 - >. D x {2u\fiW + uwi)^ 1- > "01+13 V l_ "11+13 ■ 

Ot 0*W* jg, G»M 1 go OMl,* 

We introduce the following gradings: 

|x| = -l, \t\ = -2, |«| = 1, M=3, |a|=2 

and extend them in a natural way to all polynomial functions of the internal coordinates. 
Then all computations can be restricted to homogeneous components. 



3.1 Nonlocal variables 

In subsequent computations we shall need the following nonlocal variables arising from 
conservation laws and defined by the equations 

Mx = w, 

(si)t = {-u 2 a-2uw + u 2 )/2; 
(s 3 ) x = uw, 

(s-}) t = — 2u 2 w + uu\ \ — mi«o,i; 
(se)x = w(uw-uij), 

(s 6 ), = (— u 4 a 2 — 4u i wa + 2u 2 u 2 a+4u 2 u\wi — 4u 2 w 2 —4uuo,2a + \2uu 2 w + 4uui^w 
— u\ -\-4u\u\2 — 4m i uq \w)/4\ 
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{ s i)x = —2u i W2 + 60u 2 w 2 — 36mmi.i w + 30mo,2>v + 27m 2 u\\ , 

(.s 7 ), = - I04u 4 wa -28u 4 w 2 - 132« 3 «iwi + 36u 3 u Li a - 40u 3 w 2 - 18m 2 m ,2« -48u 2 u 2 w 

— 36u 2 uiuqjCC + 144u 2 u\\w + 66u 2 uqiw\ — ISuuq^w — 36uu 2 uii + \%uu\ii\2 

— 36mki«o,iw + 30«mi.3 + 18mmj j — 36uu 2 )l a + 9wo,2«i + 18mo,2Mi,i — 30«o,3«i 
+ 36«iMo,i +3(m\u\.\uq.i — \%u\^.uqa + \Suq |W. 

The variable j,- is of grading i and computational experiment shows that for every grading i = 
4n — 2 + e, e = 0,1, there exist an such that |s,| = i. 

In addition, we found conservation laws of fractional gradings: 



(si/2)a- 


= w 1 / 2 , 


( s l/l)t 


= -w" 2 u; 


(s-ifi)* 


= w~ 3 ' 2 (6wa + w 2 ), 


( s -l/l)t 


= w~ 3 l 2 {4uwa — uw2 + 14w 2 ); 


( s -3/l)x 


= w~ 1/2 {l2w 2 a 2 + I2ww 2 a-2ww4 + 7wl), 


( s -3/2)t 


= w^^^duv^a 2 + I6uww2a + 2uww4 — 7uw2 — I24w 3 a — 32w 2 w2) 


( s -5/l)x 


= w - n/2 (216w 3 a 3 + 540w 2 w 2 a 2 - 180w 2 w 4 a + 20w 2 w 6 + 882ww2a 




— 378ww2W4 — 16ww 2 + 837iV2), 


( s -5/2)r 


= w~ n/2 (320uw 3 a 3 + 1008iw 2 w 2 a 2 - 360uw 2 w 4 a - 20uw 2 w 6 



+ 1908i(ww2a + 318uww2W4 + \6uww 2 — 837iov2 — 64u\w 2 w^,a — 1400w 4 a 2 
- 1868w 3 w2a + 580w 3 w4 + 64w 2 wiW3 -2322w 2 W2), 

etc. 

3.2 Symmetries 

A symmetry <f> = ((p w , (p u ) of Eq. d26l l must satisfy the linearized equation 

D,((p w ) +uD x {(p w ) +2ui<p K + 2wD x (<p") +wi(p u = 0, 
<p"' + D 2 (<p") - a<p" = 0. 

Direct computations lead to the following results. 

(x,t)-independent symmetries. One can observe two types of symmetries that are indepen- 
dent of x and t. The first one consists of symmetries of integer gradings: 

9" = «i; 

<P2 = uw\ +2u\w, 

<P2 = -«o,i; 

<p™ = u 2 w\0C + 2uww\ — u\w\ — 2u\\w\ — 4uQ,\wa, 
95 = u 2 ui cc + 2u 2 wi +6uu\w — u\ + 2u\2 — 2uq\w\ 
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n 

<f»g = u w\a + 2u u\wa + %u ww\ — uu\w\ + \2uu\w — 2uu\jw\ 

+ 2«o i 2M ; i — 2m[W + 2kimo j iwi +4mi j 2>v — 4mo j iw 2 , 
<jPg = 2m 3 wi +6m 2 «i>v — 3m 2 mo.iO! — 6kmo j iw — 2«o,3 +3m 2 «o,i, 

etc. Symmetries of the second type have semi-integer gradings: 

9-3/2 = (-4w 2 wia + 4w 2 W3 - 18wh'iW2 + 15wi)/(2w 7/ ' 2 ), 
<p^ /2 = -2w 1 /(w 3 / 2 ); 

(p™ 5 / 2 = (-48w 4 w 1 a 2 + 80w 4 H'3a -32w 4 w 5 -520w 3 wiW2a + 320w 3 wiW4 

+ 560w 3 w 2 w 3 + 560w 2 w] a - 1820w 2 w 2 w 3 - 2520w 2 w\w 2 2 + 6930ww 3 w 2 
-3465w?)/(12w 13 / 2 ), 

(p u _ 5/2 = {-\2w 2 w x a + 8w 2 w 3 -40ww l w 2 + 35w\)/(3w 9 l 2 ), 

etc. All symmetries are local and \(py\ = y. 

If one adds to the nonlocal setting the variables Sy (see above) then an additional series 
of nonlocal symmetries arises: 

9-1 = {s\/2W~^ 2 (4w 2 wi(X -4w 2 wt, + I8WW1W2 - 15w 3 ) + I6tw 3 w\ 

+ 2w{-4w 2 a-4ww 2 + 5wj))/(l6w 3 ), 
9-i = {si/ 2 w~ l l 2 w\ +4tu l w-2w)/(4w); 

ip™ 2 = (3s 1 / 2 w~ l/2 (-4$w 4 w l a 2 + 80w 4 w 3 a - 32w 4 w 5 - 520w 2 w l w 2 a + 320w 3 w l w 4 
+ 560w 3 w 2 w 3 + 560w 2 w] a - 1820w 2 w 2 w 3 - 2520w 2 wiw| + 6930ww]w 2 - 3465wf ) 
+ 2s _ i / 2 w~ l l 2 w i '(-4w 2 wi<x + 4w 2 wj - I8W1W2+ 15w 3 ) + 2w(96w 4 a 2 
+ 384w 3 w 2 a - 192w 3 w 4 - 14Qw 2 w\a + 1364w 2 wi w 3 +960w 2 wl 
-5514vw 2 w 2 + 3465w 4 ))/(4w 6 ), 

q>" 2 = (3s l / 2 w~ l / 2 (-l2w 2 wia + $w 2 wi - 40wwiw 2 + 35w 3 1 ) - 2s_ l / 2 w~ 1 / 2 w 3 wi 
+ 2w(24w 2 a + 24ww 2 - 35w 2 ))/w 4 , 

etc. 

(x,t)-dependent symmetries. The first three symmetries that depend on x and / are 
(pQ = t(—uw\ — 2u\w) +w, 

<p™ = t(—u 2 w\a — 2uww\ +u\w\ +2«iiwi +4uo\wa) +2(^2^1 +2uwa + u\w\), 

q>" = t(—u 2 u\a — 2u 2 w\ — 6uu\w + u\ — 2u\ t2 + 2iiQ.iw) + 2{s 2 u\ + 2uw — 2mii); 

q>4 =t(— u 3 w\a — 2u 2 u\wa — 8u 2 ww\ + uu\w\ — \2uu\w 2 + 2uu\^w\ — 2uq, 2 wi + 2u\w 

— 2u\uq\w\ — 4m 1.2 w + 41(0,1 w 2 ) +2{s 2 (uw\ +2u\w) + s^wi +4uw 2 

-4mi i iw-2m ,iwi) ) 

<p" = t(—2u 3 w\ — 6u 2 u\w + 3u 2 uq ) iCC + 6mmq,i>v + 2mo,3 — 3m 2 «o,i) +2(— s2"o,i + S3U1 
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+ u 3 a + 3u 2 w — uu 2 + 3«o,2)- 

All these symmetries, except for the first one, are nonlocal (description of the nonlocal 
variable is given in Subsection l3.1l > and, as above, the subscript denotes the grading. 

3.3 Cosymmetries 

The defining equation for cosymmetries \f/ = (ijf w , \jf") is the adjoint to the linearization 
of ED: 

D,O w ) + uD x (w w ) - m y w - y" = 0, 
2wD x (\jf w ) +wi\ff w -D 2 x {\f) + ay" = 0. 

Similar to symmetries, we consider two types of cosymmetries. 

(x,t)-independent cosymmetries. They are local and may be of integer and semi-integer 



;radings: 






= i, 


V" 


= —u\ ; 


¥4 


= u, 


¥4 


= 1*0,1', 


Wi 


= u 2 a + 2uw — u\ — 2u\ t \ , 


Wi 


= —u 2 u\ a — 2u 2 w\ — 6uu\ w + u\ — 2u\2 + 2uq^w\ 


v& 


= it 'a +4u 2 w — uu 2 — 2uu\ \ + 2«o,2 + 2w]Uo i, 


K 


= —2i?w\ — 6u 2 u\\v + 3m 2 «o,i a + 6mmq,i w + 2mq,3 ■ 



etc. and 

K /2 = W- 1 /2, 

= (4vv 2 a + 4vvw 2 -5w 2 )/(4w 3 w 1//2 ), 
\fi, 2 = 2wi/(ww 1/2 ); 
jp* = (48w 4 a 2 + 160w 3 w 2 a-64w 3 w4-280w 2 w 2 a-|-448vi-' 2 wiW3-l-336w 2 w2 

- 1848ww 2 W2 + 1 155w 4 )/(48wV /2 ), 
^ 1/2 = (12w 2 wia-8w 2 vv 3 +40wiW2-35wi)/(3w 4 w 1/2 ), 

etc. 

Similar to the case of symmetries, when one adds nonlocal variables Sy an additional 
series of nonlocal cosymmetries arises: 

Wi = {3s 1/2 w- l/2 {-4w 2 a-4ww2 + 5w 2 )-2s_ l/2 w- 1/2 w 3 +96tw 3 -10wwi)/(96w 3 ), 
yf = (-s l / 2 w^ l ^ 2 wi -4tuiw + 2w)/(4w); 
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\jr£ = (3s 1/2 w~ 1/2 (4$w 4 (X 2 + 160w 3 w 2 a - 64w 3 w 4 - 280w 2 w 2 a + 448w 2 wiw 3 
+ 336w 2 w 2 - 1848wvv 2 vv 2 + 1155w 4 ) +4.s_,/ 2 w~ 1/2 w 3 (4w 2 a + 4ww 2 -5w 2 ) 
+ I2s_i /2 w- l/2 w 6 + 2w(400w 2 w 1 a -276w 2 w 3 + 1366wwiw 2 - 1155w 3 ))/(16w 6 ), 
= (3si/ 2 w' 1 ^ 2 (12w 2 w 1 a - %w 2 W3+AQww\W2 - 35w 3 ) + 2s_ l / 2 w~ l / 2 w 3 w l 
+ 2w{-24w 2 a - 24ww 2 + 35w 2 ))/w 4 ; 
V^j = (21 s l /2 w-^ 2 (320w 6 a 3 +2240w 5 w 2 a 2 - 1792w 5 w 4 a + 512w 5 w 6 -5040w 4 w 2 a 2 
+ 16128w 4 wi w 3 a - 6912w 4 wi w 5 + \2096w 4 wja - 14592w 4 w 2 w 4 - 8832w 4 w 2 

- 8 1 3 12w 3 w 2 W2 a + 52800w 3 w 2 w4 + 181632h ,3 wiW2>V3 + 42944w 3 W2 

+ 60060w 2 w 4 a - 274560w 2 w 3 w 3 - 569712w 2 w 2 w<2 + 1021020ww 4 w 2 - 425425w?) 
+ 1 8 j_ i / 2 w~ 1 l 2 y? (48 w 4 a 2 + 1 60w 3 W2 a - 64w 3 w 4 - 280w 2 w 2 a + 448w 2 w i W3 
+ 1085_3/ 2 w _1/2 w 6 (4H' 2 a + 4ww2 -5w 2 ) +40j_5/ 2 w _1/2 w 9 + 2w(69120w 4 wia 2 
+ 336w 2 w 2 - 1 848ww 2 w> 2 + 1 1 55 w\ ) + 34 1 60w 4 w 5 + 6248 1 6w 3 wi w 2 a 

- 349880w 3 wi w 4 - 605776w 3 w 2 w 3 - 665280wVJa + 1991880w 2 w 2 W3 
+ 2772036w 2 wi w| - 7636860ww]w 2 + 3828825w^))/(72w 9 ), 

= (3i 1/ 2w _1/2 (240w 4 wia 2 -320w 4 w3a + 128w 4 >v5+2240w 3 wiW2a 

- 1344vv 3 w 1 w 4 -2240w 3 w 2 w 3 -2520w 2 w 3 a + 7728vv 2 w 2 W3 + 10416vv 2 wiW2 
-29568ww 3 w 2 + 15015wf) + 4s_ l / 2 w' 1/2 w 3 {l2w 2 w l a - 8w 2 w 3 + 4Qww\w 2 

- 35w\) + I2s_ 3/2 w' l/2 w 6 w 1 + 2w(-384w 4 a 2 - 1536w 3 w 2 a + 768w 3 w 4 

+ 3360vt.' 2 M.' 2 a-5732w 2 wiW3 -3840w 2 W2 + 23422ww 2 w 2 - 15015w 4 ))/w 7 , 

etc. 

(x,t)-dependent cosymmetries. All them are nonlocal: 
V/g v = t(-u 2 a-2uw + u 2 +2u L1 ) +2(j 2 + «i), 

y/j = t(u 2 u\U + 2u 2 w\ + 6uu\w — u\+2u\ 2 — 2«o,i w) +2(— s 2 u\ — 2uw + 2u\ t \ ); 
t//g" = f (— u 3 a — 4u 2 w + uu 2 + 2uu\ t \ — 2«o.2 — 2u\ UQ t \) + 2(s 2 u + S3 — 2«o,i ), 

= t(2u 3 w\ +6u 2 u\\v — 3u 2 uoi(x — 6mmo.iw — 2wo.3 + 3w 2 mo.i) + 2(.j 2 mo.i — ^3^1 

— w 3 a — 3w 2 w + uu 2 — 3uo :2 ) , 

etc. 

3.4 Nonlocal forms 

Recall that nonlocal forms are nonlocal variables of a special type on the ^-covering. The 
^-covering itself is obtained from Eq. d26l > by adding two additional equations 



= -tug - 2u x q w - 2wq" x - w x q u , 
= aq u -q w , 



14 



where q w and q" are new odd variables. The total derivatives on the I'-covering are 

k>Q at ik k>0 "%,k k>0 ac i\,k 



k>0 "Ik k>0 a %,k jt>0 at i\M 

The nonlocal form g; associated to a cosymmetry = (yj", Vf) (see Subsection 13. 3 b is 
defined by the equations 

A(G') = -uy"q w + (D x (wf) - 2wy?)<f - yfq\. 



3.5 Nonlocal vectors 

Dually to nonlocal forms, nonlocal vectors arise as special nonlocal variables on the £*- 
covering associated to symmetries of the initial equation. The £* -covering is the extension 
of Eq. d26b by two new equations 

p w = _ up w + UxpW+pU} 

p'', = 2wp^+w xP w + ap u , 

where p = (p w ,p u ) is a new odd variable. The total deriavatives are given by 

k>0 °Pk k>0 "P0,k k>0 °P\,k 

b, = a + I b x (-u P r + UlP » + £ P i k+l -L. + £ pt M i^r- 

k>0 "Pk k>0 "Po.k k>0 °P\.k 

The nonlocal vector P, associated to a symmetry <p = (<p w , (p u ) (see Subsection l3.2b is defined 
by the equations 

D,(P') = <p?p\ 

t> t (?) = -(uq>t +2w<p»)p w ~DM')P" ' + <PiPi- 



3.6 Recursion operators for symmetries 

The defining equations for these operators are 

D t (R w ) + ub x (R w ) + 2«iR w + 2wD x (R") + wi R" = 0, 
FT + S 2 (R") - aR" = 

(see Theorem |5), where the total derivatives are those described in Subsection 13.41 The 
following two solutions are essential: 

R*j = (g 3/2 w~ 1 1 1 ( -4w 2 wi a + 4w 2 w 3 - 1 8wwi w 2 + 15 ) 
- 4^ w 2 + lOq^wwi + q w {4w 2 a + Sww 2 - 15w 2 ))/w 3 , 
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R« , = 4(-e 3 / 2 w-'/ 2 wi +q w )/w 

and 

R w = g 3 Wl +cftwi +2q"wa, 
R» = Q 3 Ul +q w u - q\ j + g"w. 

The corresponding operators are of the form 

1 fhnw-^D- 1 ow- l / 2 /2-4w 2 D 2 x + 10w Wl D x + h & 0\ 
_I ~ 8w 3 V -4w 2 wiw- 1 / 2 Z)- 1 ow- 1 / 2 /2 + 4w 2 oj 

and 

/ W1D7 1 wiD. t + 2wa\ 
3 ~ ^mDj'+ii -Z>*+w J' 

where /112 = — 4w 2 wia + 4w 2 W3 — 18wwi>V2 + 15w 3 , h& = Aw 2 a + 8ww2 — 15w 2 . All other 
solutions obtained in our computations corresponded to operators that are generated by the 
two above. 



3.7 Symplectic structures 

Symplectic structures, as it follows from Theorem|5] are defined by the equations 

-A(S W ) -wAc(S w ) + kiS w + S" = 0, 
-2w£> x (S w ) - Wl S w + D 2 (S") - aS u = 0, 

where the total derivatives were defined in Subsection l3.4l Here are the simplest nontrivial 
solutions: 



cw _ 
5 2 — 


g 3 / 2 w~ 1/2 , 


CM _ 

5 2 — 


-q"; 


s 5 w = 


Q 3 +q"u 


s§ = 


-Q 3 ui -q w u + qi tl -q"w; 


cw _ 
5 6 — 


Q 4 + Q } u-q'{ l +q u l u-q"u u 


CM _ 

b 6 — 


-Q 4 ui + g Vl - q W u 2 - q u Q 2 + q" ( 



with the corresponding symplectic operators 

Sl -{ -l)> 

s = ( D * X D * \ 

5 \- Ul D- l -u Dx-w)' 

_ ( D^ou + uD^ 1 -D,+uD x -ui \ 

6 \-uiD~ 1 ou + UQ t \D x l — u 2 —D 2 — u 2 a — 2uw + u 2 )' 
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3.8 Hamiltonian structures 

The equations that should be satisfied by a Hamiltonian operator (see Theorem [6} are 

D,(H W ) + uD x (H w ) +2uiH w + 2wD x (H u ) + Wl H" = 0, 
H M '+52(H")-aH" = 0, 

where the total derivatives are from Subsection 13.51 In particular, we found the following 
solutions: 



hT 3 






H« 3 






i-r 2 


= 2p l w + p w\, 






= -p"; 




H" 


= P 1 w i - 2p w ww l +p"w l +2p u wa, 




H" 


= P l u\ — 2p\uw — p w {uw\ + 2uiw) - 


Pi,i+P u w. 



The corresponding Hamiltonian operators are 



1-3 



H 



-2 : 



-Dl + aD K 
D x 0J ' 

2wD x + w\ 
-1. 

W\D~ ow\ — 2ww\ w\D x + 2wa\ 

u\D~ o wi — 2uwD x — (uw\ + 2u\w) —D xt +w J 



3.9 Recursion operators for cosymmetries 

By Theorem[6] the equation to find recursion operators for cosymmetries are 

-D,(R W ) -uD x {R w ) + u x R w + R" = 0, 
-2wD x (R w ) - Wl R w +Dl{R u ) -aR" = 

with the total derivatives given in Subsection l3.5l One of solutions is presented below: 

R^=P 1 -2p w w + p" l , 

R 3 ' = -P l u\ + 2p™mv + p w (uw\ +2u\w) + p\ x -p"w 
The corresponding recursion operator is 

R 3 = 



D r owi — 2w D r 



u\D x 1 o w\ + 2uwD x + uw\ + 2u\w D xt — w 
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3.10 Interrelation 

We expose here basic facts on structural relations between the above described invariants. 
The main one is the following 

Theorem 7 Recursion operator R3 and Hamiltonian operator H_3 constitute a Poisson— 
Nijenhuis structure on the Camassa-Holm equation. Consequently, all the operators R" o 
H_3 are Hamiltonian and pair-wise compatible. In particular, R3 o H_3 = 0!H_2- 

Proof The proof consists of direct computations using the results and techniques of Ref . 1 8 0- 

A visual presentation of how symmetries are distributed over gradings is given in Ta- 
ble [2] How to prove locality of the first two series of symmetries will be discussed in Sec- 



Table 2 Distribution of symmetries over gradings 



Gradings 


-| -2 -§ -1 -i 1 2 3 4 5 6 7 8 


Local positive 
Local negative 
Nonlocal positive 
Nonlocal negative 


• • • • 

• • 

• • • • • 

• • 


tion|5] Similar presentation for cosymmetries see in Table [3] 


Table 3 Distribution of cosymmetries over gradings 


Gradings 


-§-l-|u|l§2345678 


Local positive 
Local negative 
Nonlocal positive 
Nonlocal negative 


• • • • 

• • • • 

• • • 

• • • 



The action of Hamiltonian and recursion operators for symmetries (up to a constant 
multiplier) is given in Diagram J27t : 



R_l R_ \ R_ \ R_ 1 R— 1 




The action of recursion operators for cosymmetries and simplectic structures is similar. 
We shall now prove commutativity of the local hierarchies. 



4 We are pretty sure that the pair (H_3, R-i) is also a Poisson-Nijenhuis structure and generates an infinite 
negative hierarchy of Hamiltonian operators, but could not prove this fact because of computer capacity 
limitations. 
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Lemma 1 The symmetry 93 is a positive hereditary symmetry, i.e., its action on local sym- 
metries, (p 1 — > {93,9}, coincides, up to a multiplier, with the one of the recursion opera- 
tor R3. The only symmetries that vanish under this action are <p_3/2 and q>\. In a similar 
way, the symmetry <p_i is a negative hereditary symmetry and the only symmetry that is 
taken to zero under its action is 

A direct corollary of this result is 

Theorem 8 Local positive and negative symmetries form commutative hierarchies. 



4 The scalar version 

Let us consider now the Camassa-Holm equation in its initial form (Q} with pt = and, 
similar to the matrix case introduce a new real parameter a: 

au t — u txx + 3ocuu x = 2u x u xx + uuxxx. (28) 

For the internal coordinates we choose the functions 

d'u d k+l u 
«/ = 3-7, ui k = . . , 1 = 0, 1,2, k> 1. 
dx' dx'dt" 

The total derivatives in these coordinates are of the form 

a £ a / a a 

-+Lr(u i ) — 



D x = 3- + $2«Z+13 1" E M U3 ^213 ^("3)3— 

dx £g ^ V ouQ. k dui,k ou 2 



a 2 a 2 a 



, =0 g» M/ aw/i- 

where M3 = (a«o,i —"2,1 + 3awwi — 2u\U2)/u. The equation becomes homogeneous if assign 
the following gradings: 

|*| = -1, |f| = -2, |u| = 1, \a\ = 2. 



4. 1 Nonlocal variables 

We introduce nonlocal variable associated to conservation laws of the equation at hand: 

(s 2 ) x = ua-u 2 , 

( s i)i = (—3u 2 a + 2uu2+u 2 )/2; 

i s ^)x — u{uoc — U2), 

(si) t = —2iP a + 2u 2 U2 + uu\ \ — uqiui ; 

( s 6) x = u 3 cc 2 — 2u 2 u 2 a — uu\.\OL + uu\ +M11M2, 

(j 6 ), = (— 9u 4 a 2 + I2u 3 u2(x + 6u 2 u 2 a + 4u 2 u\.ia — 4u 2 u 2 . -%uuq,\u\(x 

— 4uuq.2(X — 4mM[M2 + 4uu\U2 : i — 4uu\ : \U2 + 4uo,iu\U2 — u\ -\- 4u\u\ 2) / 4; 
( 1?7 ) r = 60m 4 a 2 — 116w 3 M2a — \2u 2 u\a — 40u 2 ui t \a-\-56u 2 u\ — 6uuQ^u\a 
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+ 28m«o,2 0! + 12m«i«2 + 10MMiK2,l +40mMi i iM2 + 2m«2,2 +2«q id! — 4«o,lMlM2 

— 2«o,lM2,l — 30mo 2"2 + 27«i«i i, 

(.57), = - 144m 5 a 2 + 240m 4 m 2 « + 48m 3 M 2 a + 124m 3 mi . 1 a - 96m 3 m| 

9 9 9 9 9 9 

— 156m «o,i«iC< — 124m mo^O! — 48m M[M 2 + 8m M1M24 — 88m u\\U2 + 2%u m 2j 2 

— 56mm ( 2 j a + 1 12m«o.i mim 2 + 38m«o.i «2,i + 78m«o,2 m 2 — 36mm 2 mi.i + \ %uu\u\^, 
+ 18mm 2 j +30m«i 3 — 18«o j i«2 + 36mo.i« 3 +36«o,i«i"i,i — 18mo,imi ^ + 9mo,2«i 

+ 18mo.2«1,1 — 30mo.3Mi. 

4.2 Symmetries 

A symmetry 9 must satisfy the linearized equation 

aD,((p) -D 2 x D,((p) - uD 3 x ((p) - 2uiD 2 x ((p) + (3au - 2u 2 )D x ((p) + (3a«i - M 3 )(p = 0. 

We computed two types of symmetries. Everywhere below the subscript was chosen in 
a way to correspond the enumeration taken for the matrix case. 

(x, t)-independent symmetries. We present the first four of them: 
(p\=u\, 

<Pl = -1*0,1, 

95 = 3m 2 miO! — 4mmo,iC< — 2uu\U2 + 2mm 2j i + 2mo i im 2 — u\ +2mi 2 , 

^999 9 

9e = 2m mi a — \\u uoia — 2u u\U2 + 2u «2,i +6m«o.i«2 +3M0.1M! — 2«o.3- 
All these symmetries are local. 

(x,t) -dependent symmetries. The first three (*,r)-dependent symmetries are 
90 = f«o,i +" ; 

9? = 2s2U\ +t(—3u 2 ui<X + 4uuo,i a + 2uu\U2 — 2uu2,\ — 2«o,i"2 + u\— 2u\2) 

+ 4(u 2 a — uu2 — «u), 
94 = — 2^2^0,1 + 2.S3M1 +f (— 2m 3 mi a + 1 1m 2 mo,i 0! + 2u 2 u\u2 — 2m 2 «2,i 

— 6mmo i iM2 — 3mo.im 2 +2mo j 3) +2(4w 3 a — 3m 2 m 2 — uu\ + 3mo,2) 

The only local symmetry in this series is 91 . 

4.3 Cosymmetries 

The defining equation for cosymmetries is 

-OCA(V) +D 2 x D,{y) +mD 3 (v/) +miD 2 (v/) + (w 2 - 3«m)/) x (9) = 0. 
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(x,t)-independent cosymmetries. These cosymmetries are local. The first four of them are 

¥3 = 1, 
Y4 = u, 

\ff-j = 3u 2 a — 1uu2 — Mi — 2u\ t i , 

\fffi = 5u 3 a —4u 2 U2 — uu\ —2uu\ \ +2mo.imi +2mo,2- 

(x,t) -dependent cosymmetries. These cosymmetries are nonlocal, except for the first one: 

\j/ 2 = —tu, 

% = 2s 2 + t(-3u 2 a + 2uu 2 + u 2 + 2u 1A ) + 2u x , 

V>6 = 2s 2 u + 2.?3 +t(— 5u 3 a + Au 2 u2 + uu\ + 2uu\ t \ — 2wo,i"i — 2^0.2) — 4mo,i, 

etc. 



4.4 Nonlocal forms 

Nonlocal forms arise in the £-covering which is given by the equation 

aq, - q xxt - uqjca - 2u x q xx + (3(Xu - 2u xx )q x + (3(Xu x - u xxx )q = 
with the total derivatives 

Ac=Ac + I>+l^- + £ (qi, k ^—+q2,k^—+D^(q 3 )^—) , 

t6 dqi V oqo,k oqi,k oqi,kJ 
2 a 2 a 

/to /totM °W 

where 93 = (ago.i — <?2.i + (3«ia — i<3)g + (3ua — 2u 2 )q\ — 2u\q2)ju. 

To any cosymmetry xj/j we associate a nonlocal form Q defined by the equations 

D x (Q i ) = (a\ l f i -D 2 x (x l f i ))q, 

A(G') = ((u 2 -3ua)\ifi + uDl(xif i ))q u + (ui\if i -uD x (xif i ))q l +uxif i q 2 
-Ac(Vi)?o,i+V'i?i,i- 

4.5 Nonlocal vectors 

Nonlocal vectors arise in the £*-covering. The latter is the extension of the initial equation 
by the equation 

-ap t + Pjoc, + up xxx + Ujcpja + (u„ - 3dll)p x = 

with the total derivatives 

A d n / d d „t. , d \ 
D x =D x +Y,Pl+i^-+Y,[ Pi-kj + d ,(pi) 



1=0 



°Pl K>\ V d PO,k OP\,k d P2,k / 
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where pi = {{3ua-u2)pi -wipi + otpo.i —p2,\)/u. 

To any symmetry <p,- there corresponds a nonlocal vector P' defined by the relations 

(P i ) x = (< Pi a-D x -(<p i ))p, 

(pi), = ((u2-3ua)<Pi + uiD x ((p i )+uDl((p i ))p-uD x ((p i )pi+u<p i p2 
-D x ((pi)p 0t i + q>ip\ t \. 

4.6 Recursion operators for symmetries 

The defining equation for these operators is 

ab,(R) - D 2 D,(R) - uD 3 x (R) - 2uiE>l(R) + (3au - 2u 2 )D x (R) + (3ce«i - m 3 )R = 

where the total derivatives are those presented in Subsection 14.41 We consider here two 
nontrivial solutions: 

R_! = (4Q 3 / 2 (ua — ii2y 5 / 2 u(2u 2 ma 2 + uuo.ia 2 — Auu\U2a — m«2,i« 

2 2 1 3 2 2 

— Mo, \U2Ol + 2mi« 2 + "2^2,1 ) — 4^2 M (w0£ — U2) + 4#i (wG! — W2) ~ u(—2u U\a 

— UUq\ OL 2 +4uil\U2(X + UU2A OL + U0/1U2OC — 2u\ U 2 — U2U2.1) 

+ 2q((ua — U2) ~ 3 {4u 2 u\a 2 + AuuQ\u\a 2 — %uu\u2<x — Auu\U2^a 

+ Mq 1 OL 2 — 4«ojMi J/2O! — 2mojM2.iO! + 4u 2 U 2 + \u\U2U2,\ +U21) 

+ 2u 2 a(ua - u 2 y 1 ))/u 2 

and 

R3 = Q 3 U\ — 171.1 — q2U — q\U\ + q(2u(X — U2). 
The operator corresponding to the second one is 

R3 = au\D~ x —D xl —uD 2 —u\D x + 2ua — ii2- 
The first operator is too complicated to present it here. 

4.7 Symplectic structures 

The equation defining symplectic structures is 

-aD,(S) +D^D t (S) + uS x (S)+uiDl(S) + (u 2 -3au)D x (S) = 0, 
where the total derivatives are from Subsection |4.4| The simplest solutions are 

s 5 = e 3 , 

S 6 = Q 4 + Q 3 u-q ,i -q\u 
while the corresponding operators have the form 

5 5 = CCD' 1 , 

5 6 = D x l o (au — 112) + auD x l —D, — uD x . 
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4.8 Hamiltonian structures 

Hamiltonian structures are defined by the equation 

aA(H) - 5?A(H) - uD 3 x (H) - 2miDJ(H) + (3au - 2u 2 )D x (H) + {3au x - m 3 )H = 0, 

where the total derivatives are given in Subsection l4.5l We present here three solutions. Two 
of them are local and the third one is nonlocal: 

H-3 = PU 

H-2 = -po.i -piu+pui, 
Hj =P u\ - pi.2 - pi,\u + po.iua - P2U0.1 + piu(-ua + u 2 ) 
+ p(—uu\a + UQ.\a + u\U2) 
The corresponding operators are 
H_3 = D x , 

H_2 = ~~ D, — uD x + «i, 
Hi = mi .D^ 1 {{—2uu\a — mo,i u + 2u\u 2 + «2,i )/ u ) — uD^D, — D x Dj + uaD, 
— uo,\D x + u(—ua + ii2)D x — Auu\a + uuj, + 3u\U2 + U2,\ ■ 

4.9 Recursion operators for cosymmetries 

These recursion operators are defined by the equation 

-aD,(R)+D 2 x D,(R) + uD 3 x (R)+ Ul D 2 x (R) + (u 2 -3au)D x (R) = 0, 
where the total derivatives are given in Subsection l4.5l One of the solutions is 

R3 = P 1 + pii + p2U — 2ua + U2 

to which the operator 

R3 = D x l o ((— 2kmi a — mo.i oc + 2u\U2 + u2,\)/u) +D xl + uD 2 — 2ua + u2 
corresponds. 

5 Conclusion 

We finish this paper with a number of remarks. 

Remark 1 First of all, let us stress again what has been done above. We treated the Camassa- 
Holm equation directly, without artificial assumptions about its evolution or pseudo-evolu- 
tion nature. The pass from the original form {TJ to system d26l l (which is not in evolution 
form) was done by technical reasons only. On this way, we found an infinite family of pair- 
wise compatible Hamiltonian structures, recursion operators for symmetries and cosymme- 
tries, and symplectic operators. These structures lead to existence of two commutative series 
of local symmetries and conservation laws and thus the equation is integrable. 

Remark 2 Several comments are worth to be made in relation to Theorem[3] Two cases must 
be distinguished: local and nonlocal. In both cases the problem reduces to reconstruction of 
nonlocal shadows up to symmetries and subsequent computation of the bracket. 
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Local case. This case is simple and the result is as follows: any shadow in the ^-covering 
corresponding to a recursion operator can be canonically lifted up to a symmetry of the 
covering equation and the commutator of these symmetries corresponds to the Nijenhuis 
bracket of the operators. A similar result is valid for the shadows in the ^-covering cor- 
responding to local bivectors (in this case the commutator of symmetries is related to the 
Schouten bracket). 

Nonlocal case. The case of nonlocal operators is much more complicated and rests on the 
problem of how to commute nonlocal shadows. Seemingly there exists no natural definition 
of such a commutator, but in Ref. |9| we proposed a procedure both to reconstruct and 
commute shadows. A weak point of this procedure is that it is based on rather cumbersome 
and intuitively not obvious notion of shadow equivalence. Dealing with equivalence classes 
of shadows necessitates elimination of certain kind of nonlocal variables that we call pseudo- 
constants and that are intrinsically related to the covering at hand. This not always simple to 
do and we plan to simplify and clarify the procedure. 

Remark 3 Of course, all invariants of the Camassa-Holm equation in its initial form can 
be obtained from those computed for the system by simple transformations. E.g., to obtain 
symmetries of Eq. d28b from those of d26b one should substitute in the M-component the 
variables w; by D l x (ua — ua). On the other way, any cosymmetry in the scalar case can 
be constructed directly from the corresponding cosymmetry of Eq. J26I I by changing the 
variables vt>/ in the w-component by D[{ua — 112). Similar transformations are applicable to 
other structures. 

Remark 4 Computation of symmetries and cosymmetries in the matrix representation can 
be simplified using the following observation: 

Proposition 1 The following correspondence between the components of symmetries and 
cosymmetries of matrix equation d26l > is valid: 

w D,+uD x -u x „ (-1) „ a-D 2 x 

Wk 1 > v4 1 ► 1 — > Ik-v 

Remark 5 Now we shall show how to prove locality of symmetry hierarchies. Let us intro- 
duce a new variable v = w 2 . Then Eq. J26l l transforms to 

v, = -u x v - uw x , 

V 2 = au-Uja, 

i.e., the first equation acquires potential fornjfl This means that for any symmetry <p = 
(<p v , <p") the form aty = <p v dx — (u<p v + v<p") dt is a conservation law of Eq. d29l l. Compar- 
ing gradings, it is easily checked that all conservation laws of the form Oty are trivial and 
consequently <p v lies in the image of D x for any symmetry (p. From this fact it follows that 
the action of recursion operators is local. 

Remark 6 Finally, we indicate relations between the Camassa-Holm equation and the equa- 
tion describing short capillary-gravity waves (the Neveu-Manna equation, see GHH)) 

u xy = u- uu xx + ^u xx u 2 x . 

We strongly believe that there exists a deformation connecting these two equations and their 
integrable structures are closely related to each other. We intend to discuss this relation 
elsewhere. 



This fact is related to existence of the nonlocal variable S1/2, see above. 
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